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ABSTRACT 

Nonlinear  systems  which  obey  a principle  of 
superposition  under  some  operation  have  been 
termed  homomorphic.  A class  of  homomorphic  systems 
which  has  found  wide  application  is  that  involving 
the  filtering  of  siqnals  that  have  been  combined  by 
convolution.  The  use  of  homomorphic  systems  to 
deconvolve  both  seismic  reflection  and  teleseismic 
data  has  been  proposed  and  explored  by  a number  of 
researchers  with  varying  success.  The  resultant 
strategies  may  be  globally  characterized  as  a 
deterministic  approach  to  signal  analysis  in  that 
no  account  is  made  for  realistic  deviations  of 
the  data  from  the  idealized  time-invariant  seismic 
models.  Several  novel  results  are  presented  in 
this  paper.  A new  class  of  systems,  called  band- 
pass homomorphic  systems,  is  discussed,  which  are 
matched  to  the  band-pass  nature  of  seismic  signals. 
The  implementation  of  homomorphic  systems  is  im- 
proved and  made  more  reliable,  through  the  use  of 
a new  phase  unwrapping  technique.  Finally,  the 
concept  of  short-time  homomorphic  analysis  is 
introduced  and  new  strategies  for  wavelet  estima- 
tion by  homomorphic  filtering  are  proposed. 


I . INTRODUCTION 

An  important  class  of  seismic  analysis  tech- 
niques is  based  on  a representation  of  the  seismic 
signal  as  a convolution  of  components,  with  the 
basic  signal  processing  task  being  deconvolution. 

There  are  a number  of  approaches  typically 
available  for  carrying  out  a deconvolution.  One  of 
the  most  common  is  the  use  of  linear  inverse  fil- 
tering; that  is,  processing  the  composite  signal 
through  a linear  filter,  whose  frequency  response 
is  the  reciprocal  of  the  Fourier  transform  of  one 
of  the  signal  components.  Obviously,  in  order  to 
use  inverse  filtering  such  components  must  be 
known  or  estimated.  Homomorphic  signal  process- 
ing [1)  represents  an  alternative  approach  to 
seismic  deconvolution.  This  class  of  filtering 
techniques  represents  a generalization  of  linear 
filtering  problems. 

The  objectives  of  this  paper  are  first  to 
review  the  various  ways  in  which  homomorphic  fil- 
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tering  has  been  explored  for  use  with  seismic  data 
processing;  and  second,  to  present  a number  of 
r.c iw  results  in  homomorphic  signal  analysis  involv- 
ing not  only  the  theory  and  implementation  of 
homomorphic  systems,  but  the  overall  analysis 
strategy  as  well.  Before  doing  so,  however,  we 
discuss  in  the  next  section  the  basic  principles 
of  homomorphic  filtering. 

II.  HOMOMORPHIC  SYSTEMS  FOR  CONVOLUTION 

Where  we  are  faced  with  the  problem  of  fil- 
tering signals  which  have  been  added,  we  very 
often  use  a linear  filter.  A principal  reason 
is  that  linear  filters  are  analytically  convenient 
for  dealing  with  signals  that  have  been  added. 

The  principal  advantage  with  linear  filtering 
applied  to  added  signals  is  that  if  the  behavior 
of  the  filter  for  each  of  the  components  separate- 
ly is  known,  then  the  behavior  for  the  sum  is  the 
sum  of  the  responses,  as  a consequence  of  the 
principle  of  superposition. 

The  notion  of  homomorphic  deconvolution  is  to 
apply  the  same  basic  reasoning  to  deconvolution. 
Towards  this  end,  we  consider  a class  of  systems 
matched  to  convolution  in  the  same  way  that  linear 
systems  are  matched  to  addition.  Then,  the  class 
of  systems  that  we  consider  have  the  superposition 
property  illustrated  in  Figure  1.  It  can  be  shown 
that  this  class  of  systems  can  be  represented  by 
the  cascade  of  three  systems  as  shown  in  Figure  2. 


Figure  1:  Convolutional  Superposition 
Property 
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Figure  2:  Canonic  Representation  of 
Homomorphic  Systems  for  Convolution 


The  first  system  , with  transformation  Dt[*]  maps  a 
convolution  of  components  to  a sum.  The  system  Ll*| 

is  a linear  system  and  the  system  D,1]”]  is  the 
inverse  of  the  system  D,I*1  and  consequently  maps 
a sum  of  components  to  a convolution.  The  system 
n4  1 • 1 and  its  inverse  are  specified  by  the  charac- 
teristics of  the  signal  space  with  which  we  are 
dealing.  Thus,  all  the  flexibility  within  this 
class  lies  in  the  linear  system  L[‘I. 

Full -Band  Homomorphic  Systems 

The  transformation  D4(-]  has  been  defined, 
more  specifically,  by  the  property  that  the  z- 
transform  (or  Fourier  transform)  of  its  output 
is  equal  to  the  complex  logarithm  of  the  z-trans- 


form 

(or  Fourier  transform)  of  its  input. 

i.e. , if: 

x (n)  * D#  lx  (n)  1 

(la) 

then 

- 

X(z)  * log  x(z) 

(lb) 

Since  the  z-transform  is  a complex-valued  function 
its  logarithm  must  be  carefully  defined.  In  par- 
ticular, since  the  transformation  D I«1  maps  con- 
volution to  addition,  the  transformation  of  eq.(lb) 
is  required  to  map  multiplication  to  addition  and 
this  imposes  a particular  interpretation  on  the 
complex  logarithm.  The  input  space  to  Dt(*]  is 
restricted  to  the  set  of  all  full-band  sequences 
x(n),  that  is,  all  stable  sequences  with  stable 
inverses,  otherwise,  the  complex  logarithm  would 
not  be  analytic  on  the  unit  circle  and  the  output 
of  D.|*]  would  be  unstable.  Then,  in  particular, 
low-pass,  high-pass  and  band-pass  filtered  sequen- 
ces are  not  allowable  inputs  to  a homomorphic 
system  with  a system  Da(*l  as  in  equation  (1).  We 
shall  refer  to  this  class  of  systems  as  full-band 
homomorphic  systems. 

The  output  xin)  of  the  system  D l*J  is  referred 
to  as  the  complex  cepstrum  of  the  input  x(n).  This 
terminology  is  motivated  by  the  relationship 
between  this  transformation  and  the  cepstrum  as 
proposed  by  Bogert,  Healy  and  Tukey  for  the  detec- 
tion of  echoes.  Specifically,  the  cepstrum  of  a 
signal  was  defined  as  the  power  spectrum  of  the 


logarithm  of  the  power  spectrum.  Since  the 
cepstrum  was  directed  toward  echo  detection  rath- 
er than  a deconvolution,  retention  of  phase  infor 
mat  ion  was  not  important.  Thus,  it  does  not  uti- 
lize phase  information  and  involves  the  logarithm 
of  real,  positive  values.  By  contrast  the  output 
of  the  system  D# [ ■ ] is  referred  to  as  the  complex 
cepstrum  in  reference  to  the  fact  that  it 
requires  the  use  of  the  complex  logarithm.  It  is 
important  to  note,  however,  that  the  complex 
cepstrum  is  a real  valued  sequence.  It  should  be 
noted  also  that  the  cepstrum  as  defined  by  Br>gert 
Healy  and  Tukey  is  proportional  to  the  even  part 
of  the  complex  cepstrum. 

Properties  of  the  Complex  Cepstrum 

In  carrying  out  a deconvolution  using  homomor 
phic  filtering,  an  appropriate  choice  for  the 
linear  system  Ll*l  must  be  made.  This  choice  is 
intimately  connected  to  a number  of  properties  of 
the  complex  cepstrum,  as  suiranarized  below. 

Property  1 - The  complex  cepstrum  of  a convo- 
lution of  two  (or  more)  signals  is  the  sum  of 
the  individual  complex  cepstra. 

Property  2 - The  complex  cepstrum  of  a 
minimum-phase  sequence  is  zero  for  n ■'  0 and 
the  complex  cepstrum  of  a maximum-phase  sequence 
is  zero  for  n > 0. 

Property  3 - The  complex  cepstrum  win)  of 
a pulse  w[n } whose  spectrum  is  smooth  tends  to 
be  concentrated  around  low-time  values. 

Let  us  denote  by  r[n]  an  impulse  train,  of 
the  form 

M 

rln)  - I rk  Sln-^l  , nk+1  > n , rk  f 0,vk 
k-1 

where  (r,  ,n,  ) k=l,...,M  denote  the  arrival  times 
k k 

and  amplitudes  associated  with  each  impulse  in 
r[nj.  Let  us  define,  in  general,  a periodic  im- 
pulse train  with  period  T > 1,  as  one  for  which 
the  interarrival  times  are  multiplies  of  T,  that 
is : 


"k  ■ Vl  + *k  • T-  k‘2 M'VJn  * ° 


then  it  can  be  shown  that: 

Property  4 - The  complex  cepstrum  r(nl  of  a 
periodic  impulse  train  r(n)  is  also  a periodic 
impulse  train  with  the  same  period,  that  is: 

CO 

r In)  - l r Ik]  6ln-kT] 
n— “ 


Property  S - Consider  a minimum-phase  im- 
pulse train  for  which  the  first  interarrival 
time  is  T , so  that 

n2  ” ni  ’ Ti'  nk  «rbitrary>  * > 2 
Then,  the  complex  cepstrum  Is  zero  for  0 < n < Tj. 


exponentially  weighted  as 


In  qeneral , r[n]  is  non-zero,  only  at  times  0, 

n -n  , n ~n  ,...,n  -n  as  well  as  at  all  positive; 

* i 1 MI 

linear  combinat ions  of  these  times.  A similar 
result  can  be  derived  for  maximum-phase  impulse 
trains. 

Property  6 - The  cepstral  structure  of  a 
mixed-phase  impulse  train  may  be  very  sensitive 
to  minor  chanqes  in  the  impulse  amplitudes.  In 
general,  the  cepstra  of  mixed  phase  impulse  trains 
exhibit  a rather  elaborate  structure,  which  with- 
in the  limits  imposed  by  our  present  understanding 
of  the  cepstral  mapping,  offer  no  clues  regarding 
the  corresponding  time  structure  and  vice  versa. 

This  characteristic  should  be  kept  in  mind,  as 
it  plays  a major  role  in  shaping  strategies  for 
data  analysis  by  homomorphic  filtering,  to  be 
described  in  Section  IV. 

III.  SEISMIC  DECONVOLUTION  BY  TIME 
INVARIANT  HOMOMORPHIC  ANALYSIS 

Seismic  signals  are  often  represented  in 
time- invar iant  form  as  the  convolution  of  a 
seismic  wavelet  w[n],  thought  to  represent  the 
combined  source  and  earth  attenuation  effects  with 
an  impulse  train  r[n),  representing  the  earth 
impulse  response.  A third  component  pin]  may  also 
have  to  be  included  to  represent  the  presence  of 
water  reverberations  or  ghosting  mechanisms. 

The  structure  of  the  complex  cepstrum  sin]  is 
such  that  the  complex  cepstrum  of  the  seismic  wave- 
let win] , that  is  the  reflector  series  ?(n]  and 
that  of  the  reverberation  train  pin], may  tend 
to  occupy  disjoint  time  intervals.  For  example,  if 
the  source  pulse  has  a relatively  smooth  spectrum, 
wfn]  will  tend  to  be  concentrated  near  the  origin. 
If  the  reflector  series  is  minimum-phase  then  it 
contributes  to  the  complex  cepstrum  only  for 
times  greater  than  or  equal  to  the  time  of  the 
first  reflection.  For  marine  seismograms,  the 
reverberation  train  is  a periodic,  minimum-phase 
impulse  series,  and  its  complex  cepstrum  will 
therefore  be  right-sided  and  periodic  with  the 
same  period. 

This  tendency  for  different  components  to 
occupy  different  time  intervals  in  the  cepstral 
domain  suggests  the  possibility  of  separating 
these  components  by  gating  (i.e.  windowing)  the 
cepstrum.  A window  around  the  time  origin,  for 
example,  would  tend  to  retain  the  contribution 
from  the  source  pulse.  This  low-time  windowing  on 
the  complex  cepstrum  corresponds  to  smoothing  (i.e. 
low-pass  filtering  the  low  spectrum). 

Cepstral  gating  was  first  explored  for  seismic 
data  by  Ulrych  (2],  and  was  directed  towards  the 
deconvolution  of  teleseismic  events.  Here,  the 
data  was  modelled  as 

sfn]  - win]  * rln]  (2) 

where  rln]  is  in  general  a mixed-phase  impulse 
train.  In  order  to  render  the  complex  cepstral 
component  associated  with  the  reflector  series  as 
simple  and  predictible  as  possible,  the  traces  were 


s' In]  = slnj  un  = (w(n]’in)  * (r(n]  »n  ) = 

w'  [n]  * r ' (nj . (1) 

prior  to  homomorphic  deconvolution.  This  has  the 
effect  of  radially  scaling  the  z-transform  of  a 
sequence.  In  particular,  any  zero  of  R(z)  outside 
the  unit -circle  in  the  z-plane  can  be  moved  inside 
the  unit  circle  by  exjxjnential  weighting,  thus 
rendering  the  weighted  impulse  train  r'(n],  mini- 
mum-phase. The  weighted  seismic  wavelet  w' [n]  was 
then  recovered  by  low-time  gating,  with  a cutoff 
frequency  which  was  determined  by  observation  of 
the  theoretical  impulse  response  of  the  earth  at 
the  recording  site,  which  was  known  beforehand. 

Ulrych' s approach  to  teleseismic  deconvolu- 
tion was  further  pursued  by  Ulrych  et.al.  (3], 
Clayton  and  Wiggins  (4),  and  Sommerville  15]. 

These  authors  were  able  to  estimate  very  reason- 
able source  pulses  for  a number  of  earthquakes 
and  nuclear ’ explosions , using  the  homomorphic 
filtering  technique.  However,  they  also  pointed 
out  a number  of  problems  encountered  in  such 
processes.  For  example,  they  recognized  the 
computational  difficulties  of  defining  the  complex 
logarithm  in  regions  of  low-spectral  amplitudes. 
They  also  questioned  the  low  time  assumption  for 
the  complex  cepstrum  of  the  source  pulse  since  in 
some  instances  the  source  pulse  had  sharp  spectral 
zeros  thus  leading  to  localized  but  sharp  varia- 
tion in  frequency. 

The  application  of  homomorphic  filtering  in 
seismic  reflection  studies  is  discussed  first  by 
Buhl,  Stoffa  and  Bryan  16]  in  the  context  of 
shallow-water  marine  dereverberation  and  source 
deconvolution.  Here  the  data  was  modelled  as: 

sin]  * w (n]  * r (n]  * p In]  (4) 

where  p[n]  denotes  a periodic  minimum-phase  rever- 
beration operator.  Exponential  weighting  was 
applied  to  make  the  reflector  series  minimum-phase. 
Letting  T^  denote  the  two-way  travel  time  within 

the  first  earth  layer  and  T denote  the  water 

w 

column  reverberation  and  assuming 

NT  < T < (N-fl)T 
w — 1 w 

it  follows  that  high-time  gating  the  complex  cep- 
strum with  a cutoff  time  T greater  than  NT  , but 

c w 

less  than  T^  will  deconvolve  the  seismic  source, 

eliminate  the  first  N multiples  in  the  reverbera- 
tion pattern  and  reduce  the  remaining  multiples  to 
at  most  1/ (N+l ) of  their  original  pulse. 

A similar  procedure  was  also  investigated  for 
deep  water  dereverberation  and  source  deconvolu- 
tion by  Buttkus  17].  Here,  of  course,  the  cepstra 
of  the  reverberation  train  will  have  high-time 
characteristics.  De reverberation  may  then  be 
accomplished  by  notch  gating  the  cepstrum  at  the 
multiples  of  T . 
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Gallemoro  [ft]  has  attempted  to  evaluate  the 
use  of  homomorphic  deconvolution  to  marine  seismic 
data  and  to  compare  its  performance  to  the  use  of 
an  optimal  tapped  delay  line  filter.  The  criteria 
for  evaluation  were  the  percent  of  multiple  ener- 
gy removed,  percent  of  reflector  distortion 
and  visual  improvement  of  the  data  when  both 
methods  were  applied  to  synthetic  seismic  data. 

The  indications  from  the  study  were  that  homomor- 
phic deconvolution  appears  to  have  greater  poten- 
tial in  shallow  water  applications  while  the 
tapped  delay  line  filter  appears  to  be  more 
efficient  for  deep  water  data.  Buttkus  [7]  has 
also  discussed  the  limitations  of  seismic  decon- 
volution by  cepstral  gating.  Using  experimental 
evidence  he  concluded  that  the  success  of  wavelet 
estimates  by  low-time  gating  war  largely  dependent 
on  the  signal  to  noise  ratio  of  the  input  data, 
the  degree  of  overlapping  of  the  cepstrum  compo- 
nents of  the  source  wavelet  and  of  the  reflectiv- 
ity function  and  the  choice  of  the  low- time  window. 


a full-band  homomorphic  system.  Thus,  it  is 
necessary  to  redefine  the  characteristic  system 
D# [ • 1 so  that  its  input  is  matched  to  the  band- 
pass characteristics  of  seismic  data.  Such 
characteristic  systems  may  be  defined  as  a cascade 
of  two  systems,  as  illustrated  in  Figure  3. 


, “I 
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He  proposed  then  to  estimate  the  wavelet  by 
analyzing  only  a windowed  segment  within  the  trace. 
The  procedure  involved  first  the  estimation  of  the 
arrival  times  by  observation  of  spikes  in  the  power 
cepstrum  (the  power  cepstrum  is  the  square  of  cep- 
strum of  Bogert,  Healy  and  Tukey)  followed  by 
notch  filtering  the  segment's  complex  cepstrum. 
Implicit  in  this  procedure  is  again  the  assumption 
that  the  reflector  series  segment  is  minimum-phase. 

Finally,  Otis  and  Smith  [9]  attempted  to 
extend  Stockham's  blind  deconvolution  technique 
[101  to  seismic  processing.  The  technique  does 
not  involve  cepstral  gating.  Here  they  considered 
a suite  of  seismograms,  with  spatially  stationary 
source  pulse  and  (minimum-phase)  spatially  non- 
stationary reflector  series.  Then  the  complex 
cepstrum  of  each  seismic  record  will  be  the  sum  of 
the  cepstrum  of  the  source  pulse,  which  is  non-var- 
iable plus  the  cepstrum  of  the  reflector  series, 
which  varies  from  record  to  record.  By  averaging 
the  complex  cepstra  of  several  records,  an  esti- 
mate of  the  non-variable  functions  can  be  obtained 
since  the  variable  functions  will  tend  to  a mean 
value.  This  procedure  has  been  referred  to  as  log 
spectral  averaging  or  cepstral  stacking.  The  es- 
timated source  wavelet  is  then  used  to  construct  a 
Wiener  inverse  filter  to  use  as  a deconvolution 
kernel.  A similar  approach  was  also  used  by 
Clayton  and  Wiggins  14}.  We  outlined  above  the 
main  issues  involved  in  the  use  of  full-band  homo- 
morphic systems  for  the  analysis  of  time-invariant 
seismic  models.  In  the  remainder  of  this  paper  we 
shall  outline  recent  results  (111  in  homomorphic 
signal  analysis  which  greatly  enhance  the  potential 
and  the  robustness  of  homomorphic  systems  in 
seismic  data  analysis. 

IV.  NEW  RESULTS  IN  HOMOMORPHIC 
SIGNAL  ANALYSIS 

IV. 1 Homomorphic  Band-Pass  Systems 

Seismic  signals  are  often  band-pass  filtered 
prior  to  sampling  and  recording,  to  discriminate 
against  additive  noise.  As  pointed  out  in  Section 
II,  band-pess  signals  are  not  allowable  Inputs  to 


Figure  3:  Characteristic  System  for 
Homomorphic  Band-Pass  Filtering 


The  first  system  BP  corresponds  to  a frequency 
scaling  operation,  that  shifts  and  stretches  the 
pass-band  of  the  signal  to  occupy  the  entire  band. 
The  second  system  is  then  the  complex  logarithmic 
system,  which  is  characteristic  for  full-band 
homomorphic  signal  analysis. 

Let  x[n)  be  a band-pass  signal,  satisfying: 

j X[ej“l  * 0 Wj  < [ w | < w2 

X[eju)J  = J (5) 

' 0 elsewhere 

and  let  x(n]  be  the  response  of  BP  to  the  input 
x[n).  Then  x[n)  satisfies: 

X(e3“)  = X (e^W) , 0 i M i " (6) 

where 

w-w 

i(ul)  =7T  - — „ < |w|  < w (7) 

2 1 1 

This  mapping  is  illustrated  in  Figure  4.  The 
system  BP  is  referred  to  as  the  band-pass  mapping 
system.  Its  role  in  homomorphic  filtering  corres- 
ponds to  the  role  of  a linear  band-pass  filter  in 
linear  filtering  in  that  all  out-of-band  informa- 
tion is  effectively  discarded. 

The  class  ot  homomorphic  systems  with  the 
characteristic  systems  of  Figure  4 is  referred  to 
as  homomorphic  band-pass  systems.  The  theory  and 
implementation  of  this  class  of  systems  can  be 
found  in  ( 11 ) . 
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Magnitude  Spectrum  of  Band-Pass  Signal  x(n) 


Frequency 

Transformation 


(b) 


<c) 

Maqnitude  Spectrum  of  Full-Band 
Signal  x[nl 


Figure  4:  Band-Pass  Mapping 


IV. 2 Phase  Unwrapping  by  Adaptive  Integration 

The  principle  computational  stap  in  using 
homomorphic  filtering  is  the  computation  of  the 
complex  cepstrum.  The  procedure  for  determination 
of  the  complex  cepstrum  generally  relies  on  an 
explicit  computation  of  the  inverse  Fourier  trans- 
forms of  the  logarithm  of  the  Fourier  transform  of 
the  band-pass  mapped  inputs.  This  is  most  easily 
carried  out  using  the  fast  Fourier  transform  algo- 
rithm to  compute  samples  of  the  Fourier  transform 
equally  spaced  in  frequency.  Since  the  Fourier 
transform  is  a comp lex -valued  function  is  logarithm 
must  be  appropriately  defined.  In  particular,  the 


phase  of  a complex  number  is  only  unique  to 
within  an  integer  multiple  of  2^  and  this  ambigu- 
ity must  be  resolved  in  computing  the  phase  of 
the  Fourier  transform.  In  order  that  the  loga- 
rithmic transformation  map  mult iplicat ion  to 
addition  as  required,  it  can  be  shown  that  the 
phase  must  be  generated  as  a continuous  function 
of  frequency.  Thus,  if  the  phase  is  computed 
modulo  2tt,  as  is  typically  the  case,  it  must  be 
"unwrapped"  to  obtain  a continuous  phase  curve. 

Simple  algorithms,  based  either  on  examina- 
tion of  the  sampled  phase  to  locate  jumps  of  2TT 
or  based  on  numerical  integration  of  the  phase 
derivative  have  proven  to  be  generally  unreal i- 
able.  As  a result,  a novel  adaptive  numerical 
integration  phase  unwrapping  technique  has  been 
proposed  [12]  which  has  proven  to  be  very  reliable. 

IV. 3 Determinism  and  Time- Invariant  Homomorphic 
Analysis 

An  important  new  result  in  homomorphic 
signal  analysis  deals  with  the  use  of  exponential 
weighting  as  a means  of  insuring  the  minimum-ph«e 
character  of  the  reflector  series.  This  practice, 
which  has  been  used  by  many  researchers  as 
discussed  in  Section  III,  is  valid  in  noiseless 
situations  but  must  be  questioned  otherwise.  As 
it  is  shown  in  (11],  by  making  the  reflector 
series  minimum-phase  and  low-time  gating  the 
cepstra,  one  is  in  effect  constraining  the  pulse 
estimate  to  be  a replica  of  the  first  arrival. 
Thus,  although  this  strategy  may  be  reasonable  in 
teleseismic  analysis,  it  is  often  not  appropriate 
in  reflection  seismology  since  the  seismic 
wavelet  characteristics  change  in  time,  due  to 
attenuation  within  the  earth.  Furthermore,  the 
data  often  contains  significant  noise  components. 

In  summary,  the  use  of  homomorphic  systems  in 
seismic  data  processing  has  been  characterized 
by  a purely  deterministic  approach  to  signal 
analysis,  in  that  no  account  is  made  for  realis- 
tic deviation  of  the  data  from  the  idealized  time- 
invariant  models.  Such  philosophy,  on  which  the 
use  of  exponential  weighting  is  based,  leads 
typically  to  results  where  the  first  arrival 
is  perfectly  resolved  as  a very  sharp  impulse — 
thus  indirectly  confirming  the  matching  of  the 
wavelet  estimate  to  the  first  arrival — the 
remainder  being  corrupted  by  noise,  with  a 
severity  that  depends  on  the  particular  time- 
varying  nature  of  the  data  and  the  signal-to-noi se 
ratio. 


IV. 4 Short-Time  Homomorphic  Analysis 

The  homomorphic  analysis  of  signals  which 
approximately  follow  a convolutional  model  on  a 
short-time  basis  is  referred  to  as  short-time 
homomorphic  analysis.  Here,  the  basic  signal 
models  are  of  the  form: 

s(n)  * [w[nj  * r[n]J  v[n]  (8) 

where  v[n]  represents  a short-time  window.  A 
first  attempt  to  model  short-time  windowing 
effects  in  the  context  of  homomorphic  signal 
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analysis  was  recently  reported  f 1 3 1 . In  general, 
we  may  model  a short-time  segment  as 

s[nl=w'[nl*r'(n]  (9) 

where  r'[n!  represents  the  windowed  reflector 
segment 

r * [ n 1 * r (nl  • v(nl  (10) 

and  w' [n]  is  defined  in  the  frequency  domain  as 

W(e]u')  - S’  (e]U')/R’  (e]“’)  (11) 

Note  that  w’[n)  depends  not  only  on  w[nl  but  also 
on  the  structure  of  r(n)  and  on  the  window  shape, 
onset  and  duration.  If  in  carrying  out  short-time 
homomorphic  analysis,  an  accurate  estimate  of  the 
seismic  pulse  w[n]  is  desired,  then  one  must 
choose  the  window  vjn)  very  carefully.  In  parti- 
cular, one  should  not  use  windows  such  that  r'(n] 
becomes  minimum-phase  (or  maximum-phase)  since 
then  w’ [nl  will  match  the  onset  (or  offset)  of  the 
segment  s[n],  which  will  not  usually  coincide  with 
a wavelet  arrival.  By  choosing  windows  which  have 
gradual  onset  and  offset  transitions,  which  are 
long  with  respect  to  the  wavelet  length  and  smooth 
with  respect  to  the  reflector  series,  one  minimizes 
the  distortion  of  the  underlying  convolutional 
model  by  the  windowing  operation.  In  that  case, 
one  may  approximately  model  s [nl  as: 

s [nl  - w[n]  * r ‘ [nl  (12) 

A number  of  short-time  windows  have  been 
investigated  [111  in  the  context  of  short-time 
homomorphic  analysis,  namely  Hamming,  Gaussian, 
Rayleigh,  Linear  and  Exponential.  It  was  found 
that  in  general,  the  low-time  cepstral  components 
of  w' [n>  and  w[nl  were  essentially  identical,  as 
long  as  the  windowed  reflector  series  r'lnl  was 
kept  mixed-phase. 

V.  A STRATEGY  FOR  SHORT-TIME 
WAVELET  ESTIMATION 

Based  on  the  above  results,  a strategy  for 
wavelet  estimation  by  short-time  homomorphic  anal- 
ysis was  devised,  which  capitalizes  on  the  sensi- 
tivity of  the  repstral  structure  of  mixed-phase 
aperiodic  impulse  trains  to  time  domain  amplitude 
perturbations.  The  procedure  consists  essentially 
in  using  different  short-time  windows  on  the  same 
seismic  segment.  Letting  vR[nl,  k*l,...,M  denote 

a set  of  M different  short-time  windows  all  defined 
on  the  same  time-interval,  and  letting  In]  denote 

the  corresponding  short-time  segments.* 

sRfn)  * [win)  * r f n] ) • vklnj  (13) 

It  follows  then  that,  if  the  windows  are  chosen  as 
described  in  Section  IV. 4: 

sR (nj  - w(n)  * rR (n)  (14) 

where  rR  (nl  denotes  the  windowed  reflector  segment: 

rR (nl  • r {nl  • vR[nl  (15) 


The  complex  cepstra  of  each  segment  is  then  of 
the  form: 

sR  (nl  = w [n)  + rR  (nl  (lf>) 

where  it  is  expected  that  the  cepstra}  structure 
of  every  impulse  train  rR|nl  is  essentially  inde- 
pendent from  every  other.  The  seismic  wavelet 
might  now  be  estimated  using  averaging  methods, 
either  in  the  cepstral  or  in  the  time  domain. 

CEPSTRAL  STACKING:  Here  the  wavelet 

estimate  w [nl  is  formed  by  averaging  all  the 
e 

complex  cepstra  sR|n): 

1 M M 

w fnj  = — £ s.  (nl  = wfnl  + — J r’(n)  (17) 

e M , k M , k 

k=l  n= 1 

This  procedure  will  give  increasingly  better 
results  as  M-*°°.  Thus  it  may  be  applied  in  a 
fully  automatic  scheme  involving  multi-channel/ 
multi-windowing  and,  within  the  limits  of  spatial 
source  stationarity,  also  multi-shot  short-time 
wavelet  estimation,  since  then  the  total  number 
of  elements  averaged  may  be  made  quite  large. 

TIME  STACKING:  Here,  each  individual  cepstra 
is  low-time  windowed,  and  mapped  back  to  the  time- 

domain  by  D^1!*),  that  is  we  compute: 

s£ln)  = d”1  [§k (n)  • 8. (n)  1 (18) 

where  £[nl  denotes  a low-time  window  in  the 
cepstral  domain.  Since,  in  general,  the  cepstrum 
of  each  impulse  train,  rR[n]  will  have  different 

contribut ions  in  the  low-time  interval,  it 
follows  that: 

sR[n)  - wtnl  * rRL[n]  (19) 

L 

where  we  defined  rR  (n J to  be: 

rRL[nl  = D<|1[r^[nl  Mn] ) (20) 

, L 

Each  low-time  component  rR  [nl  can  be  essentially 
modelled  as: 

rRL[n)  = 6 [n-m^l  + rlR(nl  (21) 

that  is,  it  consists  of  an  impulse  plus  a noise 
component,  0R[n).  We  have  observed  [111  that  the 

cepstral  sensitivity  of  a mixed-phase  impulse 
train  to  time-domain  amplitude  changes  is  mapped, 
through  the  low-time  filtering  operation,  into  a 
similar  sensitivity  regarding  the  structure  of 
the  noise  component  nR[nl.  In  other  words,  it 

seems  reasonable  to  model  the  noise  component  at 
every  instant  n as  a zero  mean  random  variable 
with  unknown  probability  distribution.  Thus,  the 

wavelet  estimate  w [n)  is  formed  by  stacking  all 
e 

the  low-time  components,  after  appropriate  time- 
domain  synchronization,  that  is: 
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given  pass-band,  may  be  neatly  handled  by  the  class 
of  homomorphic  band-pass  systems  discussed  in 
this  paper.  The  use  of  homomorphic  analysis  based 
on  time-invariant  signal  models  may  be  appropriate 
in  many  instances,  for  example,  in  shallow  water 
marine  dereverberation,  and  in  teleseismic  data 
analysis.  Otherwise,  the  time-varying  nature  of 
the  data  as  well  as  the  presence  of  in-band  noise 
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sis of  the  deconvolution  filter  by  parameteric 
linear  predictive  modelling  of  these  wavelet  com- 
ponents has  proven  to  be  a very  versatile  tech- 
nique, often  yielding  improved  results  relative  to 
Optimum  Lag  Wiener  Filtering  (11). 
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